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We propose a simple quantum mechanical equation for n particles in two dimen- 
sions, each particle carrying electric charge and magnetic flux. Such particles appear 
in (2+l)-dimensional Chern-Simons field theories as charged vortex soliton solutions, 
where the ratio of charge to flux is a constant independent of the specific solution. As 
an approximation, the charge-flux interaction is described here by the Aharonov-Bohm 
potential, and the charge-charge interaction by the Coulomb one. The equation for two 
particles, one with charge and flux (g, $/Z) and the other with {—Zq, — $) where Z 
is a pure number is studied in detail. The bound state problem is solved exactly for 
arbitrary q and when Z > 0. The scattering problem is exactly solved in parabolic 
coordinates in special cases when q^/2Trhc takes integers or half integers. In both cases 
the cross sections obtained are rather different from that for pure Coulomb scattering. 
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I. INTRODUCTION 



o 

' Field theories with Chern-Simons (CS) term in (2-|-l)-dimensional space-time admit soliton solutions carrying 

Q>^ , both electric charge and magnetic flux [1-8]. These solutions are often called CS vortices or vortex solitons, 

■ as compared with Nielsen- Olesen vortices [9] which are electrically neutral. They appear in both relativistic 
and nonrelativistic field theories, and regardless of whether the gauge field action involves both Maxwell and 

^ • CS terms or only a pure CS term. The ratio of electric charge q to magnetic fiux $ depends only on the 

I I parameters in the field theoretical model, not on the specific solution. Such solutions are not only of interest 

(— I . in field theories, but also expected to be useful in condensed matter physics. However, the interaction of these 

d ' vortex solitons is very complicated. A single soliton solution is available in analytic form only for nonrelativistic 

^ . theory and when the Maxwell term is absent. It seems difficult to find multi-soliton solutions in closed forms, 

especially when both Maxwell and CS terms are present. Therefore a simple quantum mechanical model for 
the interaction of such vortex solitons may be of interest. The purpose of the present paper is to study such a 
model. 

. _ The real CS vortices have finite sizes. The electric charge density and the magnetic fiux density (the magnetic 

■ field) depend on the specific solution. As a simple approximation, we use point-like particles to represent them 
in this paper. Both the magnetic flux and the electric charge are then confined to a region of infinitesimal 
area, in other words, to a point where the particle is located. The vector potential associated with the fiux 
is the Aharonov-Bohm (AB) potential [10]. (see also Refs. [11,12] for some more works on the subject.) 
This is responsible to the charge-fiux interaction. As for the charge-charge interaction, we make use of the 
Coulomb potential. Note that in two-dimensional space there are two kinds of Coulomb potentials. The first 
one satisfies the two-dimensional Poisson equation with point source and is proportional to Inr, where r is the 
distance between the two point charges. The second simply imitates the form of the three-dimensional one and 
is proportional to l/r. It should be remarked that the real interaction between the CS vortices may be very 
complicated, and it depends on whether the field theoretical model involves both Maxwell and CS terms or 
only a CS term. Neither of the above forms can be expected to be capable of well describing the real situation. 
Either one is in any case a rough approximation. We prefer the latter one since it is easier to obtain exact 
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solutions in this case. This is the potential adopted in the study of the so called two-dimensional hydrogen 
atom (2H) [13-18]. 

In this paper wo confine ourselves to the framework of nonrelativistic quantum mechanics. Now that the 
forms of the interaction potentials are established, we can write down an n-body Schrodinger equation for 
these particles carrying magnetic flux as well as electric charges. This is done in Sec. II. The ath particle 
has charge and flux {qa, $a)j where a = 1,2, . . . ,n. It should be emphasized that the ratio Qa/^a does not 
depend on a, as pointed out in the first paragraph. After the time variable is separated out to obtain a 
stationary Schrodinger equation, wo concentrate our attention on the two-body problem. This is separable 
into two equations. One governs the center-of-mass motion, which is free, and the other governs the relative 
motion, which is of main interest to us and is the main subject of the remaining part of this paper. It is 
remarkable that the separability of the two-body equation crucially depends on the condition qi/^i = 92/*&2- 
We then denote {qi,^i) = {q,^/Z), (52,^2) = {—Zq,—^), where Z is & nonvanishing real number. The 
relative Hamiltonian has the same form as that for a particle of reduced mass moving in the composite field 
of a vector AB potential and a scalar Coulomb one. This may be called an Aharonov-Bohm-Coulomb (ABC) 
system. Although the so called ABC system has been dealt with by numerous works [19-24] in the literature, 
the Coulomb potential considered there is a three-dimensional one. Thus the situation is quite different from 
that studied here. In other words, the model studied in the above cited works is a three-dimensional ABC 
system, while that encountered here is a two-dimensional one. 

In Sec. Ill we study the bound state problem. Bound states are possible only when Z > 0, i.e., when the 
Coulomb field represents attractive force, regardless of whether an AB potential is present. When $ = 0, the 
spectrum is just those of the 2H. The level has degeneracy 2N + 1 (iV = 0, 1, 2, . . .). If g$/27r?ic takes 
nonvanishing integers, the spectrum is roughly the same except that the ground state has energy Ei and the 
level En has degeneracy 2A'^ (A^ = 1,2, . . .) since some solutions are not acceptable. In the general case each 
level En of the 2H splits into two, each with lower degeneracy. When q^/27Thc takes half integers, however, 
some of the splitcd levels coincide and we have again a high degeneracy. The degeneracy implies that the 
system should have SU(2) symmetry in this case, as the S0(3) symmetry of the ordinary 2H [13, 16-17]. But 
this has not been explicitly proved. 

In Sec. IV we study the scattering problem. In the general case partial wave expansion in the polar 
coordinates should be employed. However, as the asymptotic form of the partial wave involves logarithmic 
distortion due to the long range nature of the Coulomb field, it is somewhat difficult to handle the partial 
wave expansion. In this paper we restrict our discussion to special cases where q^/2'jThc takes integers or half 
integers. In these cases the scattering problem can be exactly solved in parabolic coordinates, as the ordinary 
Coulomb scattering in two dimensions [25]. Note that what we use here arc parabolic coordinates on the plane, 
and thus they are quite different from the rotational parabolic coordinates used in the discussion of the ordinary 
three-dimensional Coulomb problem in the text books of quantum mechanics. The latter are also used in the 
study of the three-dimensional ABC system [19]. When <I> = the cross section is just that for the Coulomb 
scattering in two dimensions. When g$/27r?ic takes nonzero integers, the cross section gains an additional 
term which comes from the interference of the scattered wave with an additional stationary wave present in the 
scattering solution. Without the stationary wave term the solution would become meaningless at the origin. To 
the best of our knowledge, such circumstances are not encountered previously in the literature. When q^/2T:hc 
takes half integers, the result is simple but of course rather different from that for pure Coulomb scattering. 
Without the Coulomb field our results reduce to those for pure AB scattering [10-11]. The classical limit of 
the results is also discussed. 

Sec. V is devoted to a brief summary and some more remarks. 



Consider n point-like particles carrying magnetic flux as well as electric charges in two-dimensional space. 

The ath particle has mass fia, carries electric charge and magnetic flux (qa, ^a), o- = 1, 2, . . . , n. The position 
of the ath particle is denoted by Va = {xa, Ha)- As remarked in the introduction, the ratio qa/^a is independent 
of a. More precisely, we have 



We describe the charge-flux interactions among the particles by the vector AB potentials and the charge-charge 
interactions by the scalar Coulomb ones. The n-body wave function is denoted by \l/(")(t,ri, . . . ,r„). In this 
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paper we work in the domain of nonrelativistic quantum mechanics. The Schrodinger equation for the wave 
function is then 



at 



where is the Hamiltonian of the system given by 



a=l 



Va- ^Aa(ri,...,r„) 



<la(lb 



(2a) 



(2b) 



a<b 



where the second term (the Coulomb interaction) involves a double summation subject to the condition a <b, 
and A„(ri. . . . ,r„,) is the AB vector potential at the position Va- Note that all particles, except the ath one, 
contribute to A^. Thus the components of Aq are given by 



^ax(ri, • • • , Tn) 



E^b Va - Vb 
2^|ra-rJ2' 



bjta 



A„,(r„...,r„)=X:|^^-^. (2c) 

b^a ' ' 

Since the Hamiltonian ifx docs not involve t, the time dependent factor in vp'") can be separated out. Let 

(t, ri , . . . , r„) = e-*^-*/ V^"^ (*, ri , • . • , r„), (3) 
we have for ■^("^ the stationary Schrodinger equation 

i/TV*") =i;TV'^"^- (4) 

In the following we concentrate our attention on the two-body problem, since this is the only case where 
exact analysis is possible. In this case the first summation in Eq. (2b) contains two terms and the second 
contains only one. We introduce the relative position r and the center-of-mass position R defined by 

r = ri-r2, R= — , (5) 

Where M = /xi + is the total mass of the system. Note that both Ai and A2 depend only on r, it is not 
difficult to recast Ht in the form 

where r = |r|. Using Eqs. (1) and (2c), it can be shown that 

giAi = -q2A2. (7) 

Thus Eq. (6) reduces to 

where fj, = ^1^2/(111 + /X2) is the reduced mass of the system. Now Eq. (4) can be separated into two equations. 
Let 

t/'(2)(ri,r2) = Vcm(R)V'(r), (9) 

we have 
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2^Vfl^cm = i^cm^cm, (10) 



H^) = Eip, (11a) 



where 



<^>2 y . ^2 X 

Alx = --^ 2' = 2' (lie) 

and E'en! + E = Et- Eq. (10) governs the center-of-mass motion of the system, which is obviously free and will 
not be discussed any further. Eq. (11) governs the relative motion of the two particles, which is of essential 
interest to us and is the main subject of the remaining part of this paper. In the following we omit the subscript 
r of Vr- We also denote {qi, $i) = {q, ^/Z) and {q2, $2) = {—Zq, — $), where Z is a nonvanishing real number. 
The Hamiltonian (lib) can be written as 

H = -'^{v^^^Ve\-'A (12) 

where (r, 0) are polar coordinates on the xy plane, and r has been used above. 

As pointed out in the introduction, the Hamiltonian (12) is the same as that governs the motion of a charged 
particle in the combined field of a vector AB potential and a scalar Coulomb one. However, it is quite different 
from that for the so called ABC system studied in the literature, since that is a three-dimensional model while 
ours is a two-dimensional one. More precisely, in their Coulomb potential, r = (a;^ + + z^Y^'^, whereas in 
ours r = [x^ + y^Y/"^. In fact, everything is independent of z here, or, if one prefers, there is no z component 
here. 

To conclude this section we emphasize that the separability of Eq. (4) (for n = 2) crucially depends on the 
relation (7) and thus on the condition (1). 

III. BOUND STATES 

In this section we study bound states of the two-body system. These are solutions vanishing at infinity of 
Eq. (11). It is convenient to solve Eq. (11a), with the Hamiltonian written in the form of Eq. (12), in polar 
coordinates. We denote 

^=mo + ., (13) 
where toq is an integer and < < 1. Eq. (11) can be written in polar coordinates as 

We write ip as 

V'(r,e) =i^(r)e*("-™°)^ m = 0, ±1, ±2, . . . , (15) 
then R{r) satisfies the equation 



(fR IdR 



2ijlE 2iJ,Zq^ [m + v) 



i? = 0. (16) 



Now it can be shown that E > Q gives scattering solutions which will not be discussed in this section. Thus 
bound states have -E < 0. It will also become clear in the following that bound states are possible only when 



4 



Z > 0, i.e., when the Coulomb potential represents attraction. These are all familiar conclusions in the pure 
Coulomb problem in three or two dimensions. Note that the factorized form of the solution (15) itself requires 



i?(0) = 



(17) 



except for m = TOq. This is because 9 is not well defined at the origin. It will exclude some well behaved 
solutions of Eq. (16). As E <0, we introduce a dimensionless variable p defined as 



ar, a ■ 



and a new parameter 



then Eq. (16) can be written as 



n 



IdR 
dp^ p dp 



1 A {m + uf 
4 p 



R = 0. 



Now we define a new function u{p) through the relation 



R{p) = e-p/^p^"^+Mp), 



then we have for u{p) the equation 



d'^u , , , ^ du 



TO + vl 



A U = 0. 



dp2 dp 
This is the confluent hypergeometric equation. It is solved by the confluent hypergeometric function 

u{p) = CF{\m + z/| + i - A, 2|to + z/| + 1, p), 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



where C is a normalization constant to be determined below. The other solution to Eq. (22) makes i?(r) 
infinite at r = and is thus dropped. The above solution, though well behaved at r = 0, diverges when r oo; 
u{p) behaves like e'' and R{p) like e''/-^. Therefore it is still not acceptable in general. Physically acceptable 
solutions appear when E or X takes special values so that the confluent hypergeometric series terminates. This 
happens when 



|m + z/| + 2 — A ■ 



-rir 



0,1,2,..., 



(24) 



and u{p) becomes a polynomial of order tir- Prom Eq. (19) we see that this can be satisfled only when Z > 0, 
and the energy levels are given by 



2„4 



E 



pZ'q 



2n'^{nr + \m + iy\ + l/2y' 
The corresponding wave function is 

i^nrmir, 6) = C^^^e-Pl''p\^^-\F (-n„ 2|m + + 1, p) e'("o)^. 



(25) 



(26) 



There are degeneracies in the energy levels. This is why we have not attached any subscript to E. The 
degeneracy depends on the values of v and toq. The various cases are discussed as follows. 

1. 1/ = Too = 0. This is the case of a pure Coulomb problem, or the 2H. We introduce the principal quantum 
number 



(27) 



then the energy levels are written as 



5 



2„4 



E 



N 



7V = 0,1,2,.... 



(28) 



With a given N , the possible values for {nr,m) are (iV, 0), {N — 1,±1), . . ., (0, ±A^), and the degeneraey is 
(In = 2N + 1. These results are well known [13-18]. 

2. i/ = 0, mo 5^ 0. In other words, q'^/2-Khc takes nonzero integers. In this case the energy levels are roughly 
the same. However, from Eq. (26) we see that the solution with m = is not acceptable, regardless of the 
value of rir, because the radial part of the wavefunction does not satisfy Eq. (17). Therefore the ground state 
has energy Ei, and the level En has degeneracy (In = 2N {N = 1,2, . . .). 



3. 0<z^<^or^<:^<l. In this case each level of the 2H splits into two. When m > we have 



while when m < we have 



E 



N 



2n^{N + v + l/2f' 



2Ti'^{N -v + l/2f' 



N = Q,l,2,..., 



N = l,2 



(29a) 



(29b) 



The possible values of {nr,m) correspond to E^ are {N,0), {N — 1,1), . . ., (0, iV), thus the degeneracy is 
d'^ = N + 1. Those correspond to E^ are (A^ — 1, —1), {N — 2, —2), . . ., (0, —N), thus the degeneracy is 
djf = N. The difference between the case < i/ < ^ and the case ^ < i/ < 1 lies in the order of the energy 
levels. In the first case the order of the levels is 



In the second case it is 



Eq < E^ < E^ < . . . < Ej^ < Ej^ < -Ejv+i < 

E^ < E+ < E2 < ... < EJ^ < E+_^ < E]^^-^ <... 



(30a) 
(30b) 



4. = h- In other words, q^/2Trhc takes half integers. In this case we have 



^N — ^N+1 



2„4 



2r(A^ + i)' 



A^ = 0,l,2,. 



(31) 



The degeneracy of the level is + = 2N + 2. This implies that the system has higer dynamical 

symmetry than the geometrical SO (2). It is well know that the 2H possesses SO (3) symmetry, just like the 
ordinary three-dimensional hydrogen atom possesses SO (4) symmetry. It seems that the symmetry for the 
present case is SU(2), and the above energy level corresponds to the value [N + -\- |) for the Casimir 

operator of the SU(2) algebra. But this has not been explicitly proved. One can construct the Runge-Lenz 
vector in a way similar to that in the case of 2H [16]. However, the conservation of it and the closure of the 
algebra involve some difficulty due to the singularity of the AB potential at the origin. Perhaps some other 
method should be employed to deal with the problem. 

Both bound state and scattering problems of the two-dimensional Coulomb field can be solved in parabolic 
coordinates [17-18, 25]. Here we point out that the case 2. and 4. disscussed above can also be solved in 
parabolic coordinates. As no new result can be obtained, we will not discuss the solutions in detail. In the 
next section we will deal with the scattering problem. It is in these two cases that exact solutions are available. 

Finally we give the value of the normalization constant Cn^m in the wave function (26): 



n\2nr + 2\m + v\ + l)r(2|m + iy\ + 1) 



T{nr + 2\m- 



1) 



27rnr-!(2nr-|-2|m-Fz/| -|- 1) 



(32) 



IV. SCATTERING PROBLEM 



In this section we study scattering problem of the two-body system. Here the Coulomb field may be either 
attractive or repulsive. We denote k = Zq^, which may be positive or negative. For general value of mo and v, 
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one may employ the method of partial wave expansion in polar coordinates. Then the starting point may be 
Eqs. (15) and (16). However, the asymptotic form of R{r) when r — > c» involves the Inr distortion, due to the 
long-range nature of the Coulomb field. This may be more clearly seen in the following. Thus it is not easy 
to treat the partial wave expansion and to obtain the scattering cross section in a closed form. For this reason 
we confine ourselves in this paper to two special cases where exact analysis can be carried out in parabolic 
coordinates, and defer the general discussion to subsequent study. 
Consider Eqs. (11a) and (12). Let us make a transformation 

^(r, 6») = e-*('"«+'')Vo(r, 9). (33) 
The new wavefunction ipo{r,9) satisfies the Schrodinger equation with a pure Coulomb field: 

- :^V2^o - -V'o = -EV-o- (34) 

In parabolic coordinates this equation can be separated into two ordinary differential equations while Eq. (11) 
cannot be separated. The probability current density 

j = A(^*VV-VVV*) + ^^^^^^^V>Ve (35) 



can be written in terms of ipo as 



3 = ^iroVi'o-i^oVro)- (36) 



Although ipo satisfies a simpler equation, the problem does not become easier since tpo must satisfy a nontrivial 
boundary condition 

Mr,d + 27r) = e'^^''Mr,0) (37) 

such that ^l^{r,0) is single valued. Moreover, tjjo{r,6) should have proper behavior at the origin, so that ^ is 
well defined there. The latter condition also imposes a constraint on the solution. 

It is in general difficult to deal with Eq. (37). In the following we only consider two special cases. The first 
is = 0, or q^/2Trhc takes integers. In this case Eq. (37) becomes 

Mr,0 + 2Tr) = Mr, 0), (i^ = 0), (38) 

which means ipQ is single valued. This is because the first factor in Eq. (33) is also single valued in the present 
case. The second case we are to consider is = i, or q^/2nhc takes half integers. In this case Eq. (37) 
becomes 

Mr,0 + '2TT) = -Mr,0), = I)- (39) 

Though this is not convenient in polar coordinates, it may be easily treated in parabolic coordinates. 
Now we introduce the parabolic coordinates (^, r]) whose relation with {x, y) and (r, 6) are given by 

x = ^{e-v'), y^^v, (40) 

9 9 

^ = V2rcos-, r] = V2rsm-. (41) 
In these coordinates, Eqs. (38) and (39) become 

M-^,-v) = M^,v), {>^ = 0) (42) 

and 

M-^,-v) = -M^,v), {'^ = h) (43) 
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respectively, where for convenience we have used the same notation ijjo to denote the wave function in parabohc 
coordinates. It is easy to see that other values of v in Eq. (37) renders V'ol^, v) multivalued and thus are difficult 
to deal with. Though i/jq is double valued in polar coordinates in the case = ^, it becomes single valued in 
the parabolic coordinates. This is essentially because a ^77 plane covers the xy plane twice, which is obvious 
from the relation x + iy = + iry)^/2. 

In the parabolic coordinates Eq. (34) becomes 

(5| + a^)Vo + fc'(C' + ??')Vo + Wo = 0, (44) 

where 



Note that E > Q since we are considering scattering states, and /? is dimensionless. Equation (44) can be solved 
by separation of variables. Let 

Mi,r]) = v{i)w{r)), (46) 
we have for v and w the following equations: 

v" + fc^^^u + (iikv = 0, (47) 

w" + k'^rfw + I32kw = 0, (48) 

where /3i + /32 = 4/3, and primes denote differentiation with respect to argument. The general solution of Eq. 
(44) can be obtained by superposition of solutions of the form (46) over the parameter /3i. For the scattering 
problem at hand we will see, however, that a single (ii is sufficient. No superposition is necessary. Specifically, 
we are looking for solutions that have the asymptotic property 

ipo-^e'^'', for a; ^-00. (49) 

This represents particles incident in the +x direction, as is easily verified by using Eq. (36). In the parabolic 
coordinates it becomes 

^0 ~ e''=«'-''')/^ for r}^oo and all (50) 

This can be satisfied only if 

v{^) = e^^^i"" (51) 



and w(r]) has the asymptotic form 



w{j]) ^ e-'^^'^ 1'^ , for 77^00. (52) 



It is easy to verify that t'(^) given by Eq. (51) does satisfy Eq. (47) with /3i = —i. Then the constant /32 in Eq. 
(48) is given by /32 = 4/3 + i. The subsequent discussions depend on the value of v, and we should distinguish 

between the two cases = and v 

For i/ = we define a new function u{r\) by 

w{r]) = e-''''''/\{r]), (53) 

then Eq. (48) becomes 

u" - 2ikr]u' + A(3ku = 0. (54) 
On account of Eqs. (51) and (53), the condition (42) now simply means that u{ri) is an even function of r]: 

u{—r]) — u{r]). (55) 
It is easy to find the solution of Eq. (54) that satisfies this condition: 
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u{ri)=CiF{iP,^,zkr]^), (56) 

where ci is a normalization constant. Collecting Eqs. (46), (51), (53), and (56) we obtain the solution 

^0 = cie''=(«'-''')/2j^(i/3, ^,ikr]^) = cxe^'''' F{i(3, \,ikrf). (57) 

If in addition to = we have mo = 0, i.e., for a pure Coulomb potential, this is the required solution. Taking 
the limit r — > oo, and choosing the constant Ci = e^'^/'^T(l/2 — il3)/y/w, we have for tpQ the asymptotic form 

/ r-; ■a^ 1/ m , * exp(ifcr + z/31n2fcr) 

ipo ^ exp[ikx — In k(r — x)\ + jc{0) -y= , (r ^ oo) (58) 



up to the order r where 



/c(0) = r(l/2 - i/3) exp(i/?lnsin^ 9/2 - in /A) ^ ^^^^ 



The first term in the above equation represents the incident wave while the second represents the scattered one. 
Both of them are distorted by a logarithmic term in the phase due to the long-range nature of the Coulomb 
field. Despite these distortions, it can be shown that the scattering cross section is given by 

a{e) = \fcm\ (60) 
where the subscript C indicates pure Coulomb scattering. Using the mathematical formulas 

|r(±z/3)|^= J , |r(i±z/3)p = ^_, (61) 

p smh p-K cosh p-K 

we arrive at 

= i3tanh^_ 

This is the result obtained in Ref. [25]. If mp 7^ 0, i.e., if q^/2TT'hc takes nonzero integers, the solution (57) 
is in problem, however. This is because ^/'o(r = 0) = ci 7^ 0, and according to Eq. (33) V(r = 0) = cie^*™''^, 
which is not well defined since 6 is not well defined at the origin. The correct solution for mo ^ should be 

V;o = ci[e''=^F(i/3, i,ifcr/2) _ e''=''i^(i - if3, 1, -2ikr)], (63) 

where the second term in the square bracket also solves Eq. (34) with the condition (38), and does not affect 
the boundary condition (49). We have now '^o(r = 0) = and no problem arises. Due to this additional term, 
the solution now behaves at infinity like 

V'o V'in + '^l^sc + t/Jst, {r 00) (64) 

where ipi^ and V-'sc represent the incident and scattered waves which are given by the first and second terms 
in Eq. (58), respectively, and tpst represents a stationary wave which comes from the second term in Eq. (63) 
and is given by 

_ ,iSo , A2" cos(fcr + /3 In 2kr + dp - 7r/4) 



where 



argr(i-z/3). (66) 



Since the second term in Eq. (63) is in fact the s-wave term in the partial wave expansion for a pure Coulomb 
field, the logarithmic distortion in its asymptotic form mentioned before becomes clear here. Similar distortions 

appear in all partial waves regardless of whether the AB potential is present. The first term in Eq. (64) gives 
an incident current in the +x direction (when x —00). The second gives a scattered one in the radial 
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direction (the component in the direction can be ignored when r — » oo) and leads to the cross section ac{0) 
obtained above. The third term, as a stationary wave, contributes nothing to the cross section. There are, 
however, interference terms. The interference of the first term with the subsequent ones does not lead to 
physically significant results. However, the interference of the second and the third terms actually gives rise to 
an additional term in the cross section, which will be denoted by (Tx{0). The differential cross section in the 
present case is thus given by 

(n{e) = ac{e)+a^{e), (6?) 

where 

, , V /3 tanh /Jtt cos(()o + ()i — Jliisiir 0/2) 

^><^^^ = VTk \^^\ ' ^^^^ 

and 

<5i =argr(i/?). (69) 

In the neighbourhood of ^ = 0, CTx (0) oscillates rapidly and thus the total contribution in a finite (but small) 

interval of 6 may be neglected. For large 9, especially near 9 = 1:, however, CTx {0) gives considerable contri- 
bution. It is remarkable that (Jx{9) is not positive definite and thus oi{9) may become negative somewhere. 
This means that the particles move toward the origin at some directions. To the best of our knowledge, simi- 
lar results were not encountered previously in the literature. Though the differential cross section 'Ji{9) may 
become negative at some direction, it docs not cause any trouble physically because the total cross section is 
positive (actually positively infinite due to the long range nature of the potentials). Indeed, ctx {0) gives a finite 
contribution (positive or negative) to the total cross section, while crc {9) gives a positively infinite one. 
Now we turn to the case f = \- In this case we make the transformation 

w{r]) = e-'^'^'/'^'nu{'n), (70) 

then the equation for u reads 

■qu" + 2(1 - ik'n'^)u' + 2k{2P - i)rju = 0. (71) 
The condition (43) means that u{ri) is an even function of r]. The required solution can be found to be 

uiri)=C2Fiip+llikr]^), (72) 

where C2 is a normalization constant. Collecting Eqs. (46), (51), (70), and (72) we obtain the solution 

V-o = C2e*'="77i^(z/3 + i, |, ikTj^). (73) 

Here two remarks should be made. First, as a function of r and 9, ipo is double valued, so that ip is single 
valued [cf. Eq. (33) where now u = ^] . Second, as a consequence of Eq. (43) and obvious from the above 
result, we have V'o('' = 0) = here, so that is well defined at the origin. We choose 

then the asymptotic form of tpo is given by 

/ r-, ■n^ 1/ M sin0/2 exp(7A7- + /Jln2fcr) 

V'o^expzA:x-z/31nfc r-x — — ^+/ 6* — ^-^ '-, r -> oo 74 

|smt^/2| y/r 

where 

f(0\ = /?^(-^/3) exp(t/3 In sin^ 0/2 + z3^/4) 

^^^> r(l/2 + i/3) x/2fcsin2^/2 " ^ ' 

Again note that both terms are double valued. The double valueness does not cause much trouble in the 
calculation. Using the formulas (61) the cross section can be shown to be 
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= mor = 1^^- (76) 

2k sm 6/2 

This has the same angular distribution as ac{0), but the dependence on other parameters is quite different. 

If we ignore the relation k = Zq^ and treat k as an independent parameter, we may set k = in the above 
results (note that Z = is not allowed in our formalism). Then we have 

<7i(^)=0, a,{e)= \ (77) 
ZTTKsm 0/2 

These are the AB scattering cross sections for the corresponding values of v. 

Finally we point out that the cross sections (62), (67), and (76), when expressed in terms of the classical 
velocity = hk/ii instead of fc, involve h explicitly. In the classical limit, h ^ 0, f3 — n/fiVc oo (this is 
actually realized in the low energy limit), we see that ax{0) is negligible in compared with ac{0), and both 
tanh/jTT and coth/37r tend to ±1. So we have in this limit 

ac{e) = am = a2{9) = , J^^L;, ^ (78) 

2/it;^ sm 0/2 

which is the classical scattering cross section for a pure Coulomb field in two dimensions. This result implies 
that the AB potential has no significant effect in the classical limit as expected. 



V. SUMMARY AND DISCUSSIONS 



In this paper we propose an n-body Schrodinger equation for particles carrying magnetic flux as well as 
electric charges. The ratio of electric charge to magnetic flux is the same for all particles. The two-body 
problem is studied in detail. The bound state problem is exactly solved in the general case, while the scattering 
problem is exactly solved in two special cases. 

The original intention of this work is to describe the CS vortex solitons by a simple quantum mechanical 
model. If the sizes of the solitons are small, the AB potential may be a good approximation in describing the 
charge-flux interaction. On the other hand, the real charge-charge interaction may be quite complicated, thus 
the Coulomb potential used here may be questionable. If a better form V{qa,qb, |i"a — rf,|) can be found for 
the interaction potential of charge qa at Va and charge qi, at r;,, then the n-body equation may be improved 
by substituting this potential for qaqb/\i^a — ^bl m Eq- (2b). In this case the last term qiq2/r in the two-body 
relative Hamiltonian (lib) should be replaced by V{qi,q2,r). With an improved potential, the Schrodinger 
equation might become more difficult to solve, however. Therefore, the model studied in this paper, even 
though it cannot well describe the interaction of the vortex solitons, may have some interest in itself since it 
allows exact analysis to some extent. 

Several aspects of this model that need further studies may be: the dynamical symmetry of the two-body 
system, the scattering problem for general value of u, and finally, the relativistic generalization of the model. 
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